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Abstract 

We propose to measure leptonic CP phase, after neutrino mass hierarchy is established, with an upgrade of the 
PINGU detector and using atmospheric neutrino flux. The upgrade, called super-PINGU, will require a few megaton 
effective volume at 0.5-1 GeV range to distinguish 6 in the range 7r/2-37r/2 from 0 after 4 years of operation. The 
distinguishability (similar to significance) of measuring 6 depends crucially on various flux, cross-section, event re¬ 
construction (energy and angle) and flavor identification uncertainties. We explore effects of these uncertainties on 
the distinguishability of measuring CP phase and suggest possible ways to minimize their impact. 
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1. Introduction 

Establishing CP violation in the leptonic sector is an 
outstanding problem in particle physics. Atmospheric 
neutrino flux measurements in large water/ice detectors 
can be used to determine the Dirac CP phase 6. Informa¬ 
tion on different S values is encoded in the neutrino os¬ 
cillation probabilities after propagation inside the earth 
(matter effect). A systematic shift of the probabilities 
with increasing d in a wide, ~ 0.2-2 GeV, energy range 
is key to measure S with atmospheric neutrino flux (T] . 

It has been found recently that the PINGU and ORCA 
detector with ~ 3 GeV threshold will have good sensi¬ 
tivity to determining the neutrino mass hierarchy O O 
mm. However, measurement of S would require a de¬ 
tector with larger effective volume and improved char¬ 
acteristics in the < 2 GeV range and in this context a 
future upgrade of PINGU (and also of ORCA), called 
Super-PINGU, was proposed in Ref. IH with detailed 
estimate of sensitivity. Here we present highlights from 
that work. 
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2. Methodology and Results 

To measure CP phase with atmospheric neutrinos, we 
calculate neutrino events in Super-PINGU by varying 6 
and compare with d = 0, keeping all other oscillation 
parameters fixed. We use an effective mass (both for 
and Ve) parameterized as 

pV,«{Ey} = 2.6 [log(£,/GeV) + Mt, (1) 

which can be realized for a total of 126 strings and 60 
DOMs per string. This gives an effective mass of ~ 2.8 
Mt at ~ 1-2 GeV, which is 4 times larger than PINGU. 
The number of neutrino events for a particular flavor 
a = e,N with energies and zenith angles in small bins 
A(Ey) and A(cos 0^) marked by subscript ij can be cal¬ 
culated as 

Nij^a = 2nNApT j dcosO^ 

JAi cos 6z 

f dEy Ve«AEv)daiEy, ey. (2) 

JAjEy 

Here T is the exposure time, Na is the Avogadro’s num¬ 
ber. The density of events of type a, da (the number 
of events per unit time per target nucleon), is given by 
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Figure 1: Distributions of Sfj for events in case of different 6 

values and d = 0 in the Ey-cos 6^ plane. Ey and 6^ are reconstructed 
neutrino energy and direction after smearing of the true energy and di¬ 
rection using PINGU reconstruction functions with a factor ^/3 better 
resolutions. Normal mass hierarchy is assumed. 


da(Ey,6z) = CTa^a + in terms of fluxes at the 

detector, and Oq, = + ^^Pea with correspond¬ 

ing oscillation probabilities and Pea • The original 
muon and electron neutrino fluxes at the production are 
0)0 = ey and O® = (D 0 (£„ 0,). 

We have computed the distributions of -f and 
Ve + Ve events for d = 0 and d 9 ^ 0 and take difference 
between the distributions in the Ey-cos 6^ plane to study 
their properties. Since there are errors associated with 
reconstructing the true neutrino energy and directions, 
we smear the ideal distributions with the energy and an¬ 
gular resolution functions of the detector to mimic the 
real situation. To estimate the sensitivity of measuring a 
CP phase different from zero we employ a distinguisha- 
bility parameter defined as 

Sij(f) = (Nfj-Nfj)l(rij(f), (3) 

where Nf. and are the reconstructed number of 
events in the i j-th bin in the Ey-cos 6^ plane for 6 and 
d = 0 , respectively, and cr^j(f) = N^j-\-(fN^j)^ is the total 
error in the ij-th bin. Parameter / is a measure of un¬ 
correlated systematic errors lO and we take / = 2.5%. 


Figure 2: The same as Fig. [T]but for Ve + Ve events. 


The total distinguishability 



(4) 


is a quick measure of significance. 

Figure shows the distinguishability S ij for -f 
events with / = 0 for different S values and using 1 - 
year of Super-PINGU data. Normal mass hierarchy is 
assumed. The shape of the distributions, specially their 
domain structures, is largely explained as due to grids 
of solar, atmospheric and interference magic lines in the 
Ey-cos plane. The oscillation probabilities are inde¬ 
pendent of S along these lines, thus separating regions 
of same sign distinguishability. Figure]^ shows Stj dis¬ 
tributions for Ve + Ve events. 

The uncertainties associated with atmospheric neu¬ 
trino flux, vN cross section, effective volume, etc. affect 
neutrino event distributions across bins in the Ey-cos 6^ 
plane. We include effects of these correlated uncertain¬ 
ties in our calculation with analogy to the pull method 
in;^^ analysis. In particular we minimize the following 
distinguishability parameter 
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where are the pull variables and are their standard 
values. The event distributions with varying are cal¬ 
culated as 


= aziiEjlG^yr 

x[l+AO.5 + cos0,)W^/^f), (6) 

where a is the overall normalization factor with the er¬ 
ror (Ta = 0.2, zi is the flux (flavor) ratio uncertainty 
{Ze = 1 for Ve events), with the error = 0.05; r] is the 
energy tilt parameter with (T;^ = 0.1; yS is the zenith angle 
tilt with cTp = 0.04. Figure!^ shows the minimized 
over (^k) for different correlated uncertainties as well as 
for no correlated uncertainties. A threshold energy of 
0.5 GeV has been assumed. Note that the contributions 
of Ve and channels to S are comparable. 


3. Discussion 

We estimate that after 4 years of operation and 2.5% 
systematics, Super-PINGU with 0.5 GeV threshold will 
be able to distinguish S = 7r/4,7r/2,7r, 37r/2 from zero 
with (7r/4) = (1 - 3), (7r/2) = (3 - 8), (tt) = 

(6 - 14), S^^\3n/2) = (3 - 8). The ranges depend on 
effects of different correlated systematics. These val¬ 
ues are a factor 4-6 improvement over the sensitivity of 
PINGU to S with 3 GeV threshold. 

The sensitivity of Super-PINGU to S can be further 
improved with following possibilities: 

• Decrease of energy threshold to 0.2 GeV from 0.5 
GeV with a denser array. This may increase sensi¬ 
tivity by 30%. 

• Stringent kinematical cut can be used to create a 
high-quality event sample with better reconstruc¬ 
tion of the neutrino energy, direction and flavor. 

• An increased exposure time will also increase the 
sensitivity to CP by a factor cx 

• Improved flavor identification at low energies. 

• Increase the density of DOMs or photocathode 
coverage. 



Figure 3: Distinguishability of measuring CP phase in the and Vg 
channels, with effects of various correlated uncertainties. 2.5% uncor¬ 
related uncertainty has been assumed in all cases. 
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• Statistical separation between the neutrino and an¬ 
tineutrino events. 


Our results show that Super-PINGU can be compet¬ 
itive to other proposals for measuring the leptonic CP 
phase associated with long base-line (LBL) accelerator 
experiments. 










